home *** CD-ROM | disk | FTP | other *** search
/ NeXTSTEP 3.3 (Developer)…68k, x86, SPARC, PA-RISC] / NeXTSTEP 3.3 Dev Intel.iso / NextDeveloper / Headers / ansi / m68k / math.h < prev    next >
C/C++ Source or Header  |  1993-01-07  |  5KB  |  159 lines
  1. /* Copyright (c) 1992 NeXT Computer, Inc.  All rights reserved.
  2.  *
  3.  *    File:    <ansi/m68k/math.h>
  4.  */
  5.  
  6. /* Inline routines. */
  7.  
  8. #ifndef _ANSI_M68K_MATH_H
  9. #define _ANSI_M68K_MATH_H
  10.  
  11. #ifdef INLINE_MATH
  12.  
  13. #include <errno.h>
  14.  
  15. #define acos(x) _builtin_acos(x)
  16. static inline double __const__ acos(double x) {
  17.     double v;
  18.     asm("facosx %1,%0": "=f" (v): "f" (x)); if (v!=v) errno=EDOM; return v;
  19. }
  20. #define asin(x) _builtin_asin(x)
  21. static inline double __const__ asin(double x) {
  22.     double v;
  23.     asm("fasinx %1,%0": "=f" (v): "f" (x)); if (v!=v) errno=EDOM; return v;
  24. }
  25. #define atan(x) _builtin_atan(x)
  26. static inline double __const__ atan(double x) {
  27.     double v; asm("fatanx %1,%0": "=f" (v): "f" (x)); return v;
  28. }
  29. #define cos(x) _builtin_cos(x)
  30. static inline double __const__ cos(double x) {
  31.     double v; asm("fcosx %1,%0": "=f" (v): "f" (x)); return v;
  32. }
  33. #define sin(x) _builtin_sin(x)
  34. static inline double __const__ sin(double x) {
  35.     double v; asm("fsinx %1,%0": "=f" (v): "f" (x)); return v;
  36. }
  37. #define tan(x) _builtin_tan(x)
  38. static inline double __const__ tan(double x) {
  39.     double v; asm("ftanx %1,%0": "=f" (v): "f" (x)); return v;
  40. }
  41. #define cosh(x) _builtin_cosh(x)
  42. static inline double __const__ cosh(double x) {
  43.     double v;
  44.     struct { unsigned pad:19; unsigned ovfl:1; unsigned pad2:12;} fpsr;
  45.     asm("fcoshx %1,%0": "=f" (v): "f" (x));
  46.     asm("fmovel fps,%0": "=g" (fpsr)); if (fpsr.ovfl) errno = ERANGE;
  47.     return v;
  48. }
  49. #define sinh(x) _builtin_sinh(x)
  50. static inline double __const__ sinh(double x) {
  51.     double v;
  52.     struct { unsigned pad:19; unsigned ovfl:1; unsigned pad2:12;} fpsr;
  53.     asm("fsinhx %1,%0": "=f" (v): "f" (x));
  54.     asm("fmovel fps,%0": "=g" (fpsr)); if (fpsr.ovfl) errno = ERANGE;
  55.     return v;
  56. }
  57. #define tanh(x) _builtin_tanh(x)
  58. static inline double __const__ tanh(double x) {
  59.     double v; asm("ftanhx %1,%0": "=f" (v): "f" (x)); return v;
  60. }
  61. #define exp(x) _builtin_exp(x)
  62. static inline double __const__ exp(double x) {
  63.     double v;
  64.     struct { unsigned pad:19; unsigned ovfl:1; unsigned pad2:12;} fpsr;
  65.     asm("fetoxx %1,%0": "=f" (v): "f" (x));
  66.     asm("fmovel fps,%0": "=g" (fpsr)); if (fpsr.ovfl) errno = ERANGE;
  67.     return v;
  68. }
  69. #define log(x) _builtin_log(x)
  70. static inline double __const__ log(double x) {
  71.     double v;
  72.     if (x<0.0) errno=EDOM; asm("flognx %1,%0": "=f" (v): "f" (x)); return v;
  73. }
  74. #define log10(x) _builtin_log10(x)
  75. static inline double __const__ log10(double x) {
  76.     double v;
  77.     if (x<0.0) errno=EDOM; asm("flog10x %1,%0": "=f" (v): "f" (x)); return v;
  78. }
  79. #define sqrt(x) _builtin_sqrt(x)
  80. static inline double __const__ sqrt(double x) {
  81.     double v;
  82.     if (x<0.0) errno=EDOM; asm("fsqrtx %1,%0": "=f" (v): "f" (x)); return v;
  83. }
  84. #define fabs(x) _builtin_fabs(x)
  85. static inline double __const__ fabs(double x) {
  86.     double v; asm("fabsx %1,%0": "=f" (v): "f" (x)); return v;
  87. }
  88. #define atanh(x) _builtin_atanh(x)
  89. static inline double __const__ atanh(double x) {
  90.     double v; asm("fatanhx %1,%0": "=f" (v): "f" (x)); return v;
  91. }
  92. #define expm1(x) _builtin_expm1(x)
  93. static inline double __const__ expm1(double x) {
  94.     double v; asm("fetoxm1x %1,%0": "=f" (v): "f" (x)); return v;
  95. }
  96. #define log1p(x) _builtin_log1p(x)
  97. static inline double __const__ log1p(double x) {
  98.     double v; asm("flognp1x %1,%0": "=f" (v): "f" (x)); return v;
  99. }
  100. #define rint(x) _builtin_rint(x)
  101. static inline double rint(double x) {
  102.     double v; asm("fintx %1,%0": "=f" (v): "f" (x)); return v;
  103. }
  104. #define logb(x) _builtin_logb(x)
  105. static inline double logb(double x) {
  106.     double v; asm("fgetexpx %1,%0": "=f" (v): "f" (x)); return v;
  107. }
  108.  
  109. #endif    /* INLINE_MATH */
  110.  
  111. /* Constants. */
  112.  
  113. #if    defined(__STRICT_ANSI__) || defined(_NEXT_SOURCE)
  114. #define HUGE_VAL (1e999)
  115. #endif  /* __STRICT_ANSI__ || _NEXT_SOURCE */
  116.  
  117. #ifndef __STRICT_ANSI__
  118. #define HUGE    HUGE_VAL
  119. #define NAN    (HUGE/HUGE)
  120.  
  121. #define MAXCHAR ((char)0x7f)
  122. #define MAXSHORT ((short)0x7fff)
  123. #define MAXINT    ((int)0x7fffffff)    /* max pos 32-bit int */
  124. #define MAXLONG ((long)0x7fffffff)
  125.  
  126. #define MINCHAR ((char)0x80)
  127. #define MINSHORT ((short)0x8000)
  128. #define MININT     ((int)0x80000000)    /* max negative 32-bit integer */
  129. #define MINLONG ((long)0x80000000)
  130.  
  131. #define MAXFLOAT ((float)3.4028234663852886e38)
  132. #define MINFLOAT ((float)1.4012984643248171e-45)
  133. #define MAXDOUBLE ((double)1.7976931348623157e308)
  134. #define MINDOUBLE ((double)4.9406564584124654e-324)
  135.  
  136. #define LN_MAXFLOAT ((float)8.872283935546875e1)
  137. #define LN_MINFLOAT ((float)-1.032789306640625e2)
  138. #define LN_MAXDOUBLE ((double)7.0978271289338397e2)
  139. #define LN_MINDOUBLE ((double)-7.4444007192138122e2)
  140.  
  141. #define M_E    2.7182818284590452354
  142. #define M_LOG2E    1.4426950408889634074
  143. #define M_LOG10E    0.43429448190325182765
  144. #define M_LN2    0.69314718055994530942
  145. #define M_LN10    2.30258509299404568402
  146. #define M_PI    3.14159265358979323846
  147. #define M_PI_2    1.57079632679489661923
  148. #define M_PI_4    0.78539816339744830962
  149. #define M_1_PI    0.31830988618379067154
  150. #define M_2_PI    0.63661977236758134308
  151. #define M_2_SQRTPI    1.12837916709551257390
  152. #define M_SQRT2    1.41421356237309504880
  153. #define M_SQRT1_2    0.70710678118654752440
  154. #endif /* __STRICT_ANSI__ */
  155.  
  156. #endif /* _ANSI_M68K_MATH_H */
  157.  
  158.  
  159.